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Abstract. We study the potential of a future, large-volume photometric survey to con¬ 
strain the axion mass m a in the hot dark matter limit. Future surveys such as Euclid will 
have significantly more constraining power than current observations for hot dark matter. 
Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter 
than ~ 0.15 eV decouple before the QCD epoch, assumed here to occur at a temperature 
Tqcd ~ 170 MeV; this leaves an axion population of such low density that its late-time 
cosmological impact is negligible. For larger axion masses, m a > 0.15 eV, where axions re¬ 
main in equilibrium until after the QCD phase transition, we find that a EUCLID-Iike survey 
combined with Planck CMB data can detect m a at very high significance. Our conclusions 
are robust against assumptions about prior knowledge of the neutrino mass. Given that the 
proposed IAXO solar axion search is sensitive to m a < 0.2 eV, the axion mass range probed 
by cosmology is nicely complementary. 
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1 Introduction 

Interest in the low-energy frontier of particle physics has surged over the past few years, 
especially in the area of axions and axion-like particles (ALPs) [1—4]. The ADMX microwave 
cavity experiment is getting ready for its final push towards covering a well-motivated range 
of axion masses at the sensitivity required to find these elusive particles if they were the dark 
matter in our galaxy [5-9]. New ideas to broaden the range of search masses are vigorously 
discussed [10-20]. In Korea, an entire Center of the Institute for Basic Science, the Center for 
Axion and Precision Physics (CAPP), has been dedicated to the search for axion dark matter 
and related topics in precision physics [21]. At DESY, first steps have been taken towards 
the largest ever photon-regeneration experiment to search for ALPs [22]. At CERN, the 
CAST experiment to search for solar axions has almost completed its original mission [23- 
28], and has paved the way for proposing a much larger next-generation axion helioscope, the 
International Axion Observatory (IAXO) [29, 30]. In no small amount, our present study is 
motivated by the IAXO proposal. 

Axion helioscopes such as CAST and IAXO take advantage of the two-photo axion 
interaction vertex both as a source of axions in the sun via the Primakoff effect and for the 
back-conversion of axions into X-rays in a dipole magnet oriented towards the sun. However, 
this back-conversion is not effective if the photon-axion mass difference is too large. Therefore, 
in the {m a , 5 a 7 )-P ararn eter plane, where g ar/ is the two-photon axion coupling strength, it has 
been traditionally difficult to probe the “axion line” on which the QCD axion models— 
as opposed to generic axion-like particles (ALPs)—reside. CAST has recently succeeded 
in touching the axion line in a narrow m a range around 1 eV [23]; it is the purpose of 
the proposed IAXO to push to greater sensitivity and in this way explore realistic axion 
parameter space. In this connection, we note that, independently of model dependence of 
the relationship between g nr , and m a , cosmological hot dark matter bounds can provide a 
guide as to the largest axion masses that are worth exploring with such an experiment. 
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Current observations of the cosmic microwave background (CMB) anisotropies and the 
large-scale matter distribution already impose severe constraints on the hot dark matter 
abundance; in the axion interpretation, this constraint corresponds to an axion mass limit 
of around 1 eV [31-36]. In the future, large-volume surveys such as the Euclid mission [37] 
will have a significantly enhanced sensitivity to hot dark matter [38, 39]. The aim of the 
present work, therefore, is to examine just what such a survey will be able to do for axion 
physics. 

To achieve this goal, we begin in section 2 with a brief review of axion cosmology. In 
section 3 we set up the cosmological framework and in particular discuss the used data sets 
and the cosmological model space. The actual numerical analysis for different assumptions 
about the hot dark matter contribution of neutrinos is provided in section 4. We summarise 
and discuss our findings in section 5. 


2 Axion cosmology 


2.1 Hot and cold axion populations 


The primary parameter characterising invisible axion models is the axion decay constant f a . 
or, equivalently, the energy scale at which the Peccei-Quinn symmetry is spontaneously 
broken. Mixing with the mesons at low energies induces a mass for the axion, given 

approximately by 


z 1 ' 2 U _ 0.58 eV 

~~ TTz Ya ™' K ~ fa/10 7 GeV 


( 2 . 1 ) 


where z = m u /md is the up/down quark mass ratio. The value of z lies in the range 
0.38-0.58 [40], which leads to an uncertainty in the relationship between m a and f a of less 
than ±5%. Therefore, we may equally use m a as a characteristic parameter. 

Figure 1 shows the “Lee-Weinberg curve for axions”, pieced together from several dif¬ 
ferent production mechanisms [41]. Non-thermal production by the re-alignment mechanism 
and possibly by the decay of cosmic strings (CS) and domain walls (DW) produces a non- 
relativistic population that can account for the cosmic cold dark matter. The thick black 
line in figure 1 shows the present-day cold axion energy density and its nominal uncertainty, 
assuming the Peccei-Quinn symmetry is broken after inflation so that each causal patch at 
the QCD epoch has a unique initial misalignment angle 0; when the axion field begins to 
oscillate. We follow reference [42] concerning the contribution from strings and domain-wall 
decay; other authors find a somewhat smaller contribution [41, 43]. 

If the Peccei-Quinn symmetry is broken before inflation such that 0, is the same in the 
entire observable universe, the axion abundance depends on 0? and in this sense is a random 
number [41]. Figure 1 shows some examples for the cosmic axion abundance in this scenario as 
thin black lines for the indicated 0i values. Note also that because the axion field in this case 
is present and massless during inflation, the quantum fluctuations it picks up independently 
of the inflaton field can later show up as isocurvature fluctuations in the CMB anisotropies 
unless the energy scale of inflation is sufficiently small [44-47]. This scenario would have been 
excluded if the tensor-to-scalar ratio r had been as high as initially claimed by the BICEP2 
team [48], because r = 0.2 corresponds to an unacceptably high inflation energy scale of 
~ 10 14 GeV [49-51]. However, with the results of the recent combined analysis of BICEP2, 
Keck Array and Planck data [52] being consistent with r = 0, reports of its demise may have 
been a bit premature. 
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Figure 1. Present-day axion dark matter density as a function of m a . The thermal axion population, 
which forms hot dark matter, is represented by the thick blue line. The thick black line denotes the cold 
axion population in the scenario in which Peccei-Quinn symmetry breaking occurs after inflation so 
that the visible universe contains many patches of different initial axion-field misalignment angles 0;; 
the energy density shown here subsumes both contributions from the re-alignment mechanism and 
from cosmic string (CS) and domain-wall (DW) decay according to reference [42]. The thin black 
lines pertain to the case in which Peccei-Quinn symmetry breaking occurs before inflation so that one 
single initial misalignment angle @i pervades the entire visible universe; the cold axion populations for 
several different values of 0, as indicated at the lines are shown. (Figure adapted from one supplied 
by Javier Redondo.) 


In addition, thermal processes produce axions which today play the role of hot dark 
matter [53]. The two populations co-exist because thermalisation of the cold population is 
extremely slow thanks to the axion derivative interactions with ordinary matter and photons, 
whose rates are suppressed at low energies. The thick blue line in figure 1 shows the hot 
dark matter density arising from thermal axion production. For axion masses m a > 0.15 eV, 
axion-pion interaction is the main production mechanism [54], and freeze-out occurs af¬ 
ter the QCD epoch. However, as rn a decreases and the freeze-out temperature approaches 
Tqcd ~ 170 MeV, the freeze-out epoch suddenly jumps to a much higher temperature because 
axion interactions with gluons and quarks before confinement are much less efficient [55-57] . 
Furthermore, any thermal axion population frozen out before the QCD epoch is necessarily 
diluted by entropy production during the QCD epoch, which causes a sharp drop in the axion 
energy density in the transition region at m a ~ 0.15 eV, as shown in figure 1. The exact 
transition between the pre- and post-QCD freeze-out regimes is not known, so we leave a gap 
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in the blue curve. For very large axion masses, m a > 20 eV, the axion decay a —> 2y is fast 
relative to the age of the universe, so that the hot axion population disappears. 

The contribution of thermally produced axions to the energy budget of the universe and 
hence the axion mass m a can be constrained by observations of CMB anisotropies and the 
large-scale structure distribution in the same way as we constrain neutrino hot dark matter 
and hence the neutrino mass sum. Our group has previously published limits on m a in a 
series of papers based on a sequence of cosmological data releases [31-35]; our most recent 
limit is m a < 0.67 eV at 95% C.L. for a minimal cosmological model and using Planck-era 
cosmological data [35], a number comparable to the results of other authors [36]. Hot dark 
matter constraints do not apply to axion masses exceeding ~ 20 eV because the aforemen¬ 
tioned a —> 2y decay removes the axion population. However, we note that axions masses 
up to ~ 300 keV are nonetheless cosmologically forbidden owing primarily to modifications 
to the primordial deuterium abundance instigated by the decay photons [58]. 

2.2 Axion decoupling temperature 

The freeze-out history of the axion-pion interaction at temperatures below the QCD phase 
transition was studied in detail in reference [31]. This treatment is technically valid for 
?i7. a <1 0.15 eV and therefore adequate relative to the sensitivity of the present generation 
of cosmological probes because the current upper bound is m a < 0.67 eV [35]. However, an 
extension to lower axion masses would be useful in anticipation of the next generation of 
cosmological probes. This is especially so in view of the forecasted sensitivity of the Euclid 
mission to neutrino masses, a(m u ) ~ 0.01 eV [39]; considering the similar phenomenology, a 
decoupling model for axions down to similar masses should be contemplated. 

Axions with very small masses decouple at temperatures above the QCD phase tran¬ 
sition, where the axion interaction with free quarks and gluons can be studied perturba- 
tively [56]. Perturbative treatments fail, however, when decoupling occurs close to the QCD 
epoch, and no full non-perturbative calculation exists to date. Therefore, to model axion 
decoupling in the mass range 0.01 eV < m a < 1 eV, we use the following prescription. 

1. We fix the QCD phase transition temperature Tqcd at a specific value. A typical choice 
might be Tqcd = 170 MeV, which conforms with the most recent calculations from 
lattice gauge theory [59, 60]. The chosen value of Tqcd does have an impact on the 
axion mass range that can be probed using large-scale structure observation, and we 
shall return to this point later. 

2. For large axion masses ( m a > 0.2 eV) we continue to use the same decoupling calcula¬ 
tion as in previous works. We extend this calculation down to an axion mass at which it 
returns a decoupling temperature matching Tqcd- For the choice of Tqcd = 170 MeV, 
this corresponds to m a = 0.145 eV. For all axion masses smaller than this limit, freeze- 
out occurs at T > Tqcd- 

3. For all axion masses smaller than this limiting value, we assume an effective number 
of entropy degrees of freedom of g* = 80 at the time of freeze-out, independently of 
the exact axion mass. This approximation is reasonable because g*{T) varies only 
logarithmically with temperature from g * ~ 60 at Tqcd ~ 170 MeV to g* ~ 100 at the 
electroweak phase transition Tew ~ 100 GeV, in contrast with the step-like variation in 
the spate of a few MeVs during the QCD phase transition (see figure 3 of reference [61]). 
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Figure 2. Present-day axion number density as a function of the axion mass (upper horizontal axis) 
and of the axion decay constant (lower horizontal axis), assuming Tqcd = 170 MeV (solid black line), 
Tqcd = 150 MeV (dashed blue line) and Tqcd = 190 MeV (dash-dotted red line). 


As we shall see later, axions decoupling at temperatures above Tqcd are practically invisible 
to the surveys discussed here, whereas axions decoupling after the QCD epoch have a clear 
impact on cosmological observations. 

Figure 2 shows the present-day axion number density as a function of the axion mass 
according to the model described above. For the choice of Tqcd = 170 MeV, the number 
density plummets from n a = 36.45 cm -3 to 10.04 cm -3 at m a = 0.145 eV, reflecting the 
jump of the axion decoupling temperature when the universe transits through the QCD 
epoch. Since the dividing line hinges strongly on the assumed value of Tqcd 5 we show 
also the cases of Tqcd = 150 and 190 MeV, where the corresponding drops in n a occur at 
m a = 0.185 and 0.117 eV, respectively. Given the uncertainty in Tqcd and the lack of a full 
quantitative calculation of the axion freeze-out in this temperature region, we cannot state 
the sensitivity to axion masses any more precisely. Henceforth, we shall use Tqcd = 170 MeV 
as our benchmark case, and 

tbqcd = 0.145 eV (2.2) 

shall denote the limiting axion mass corresponding to nominal freeze-out at this temperature. 
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3 Cosmological set-up 
3.1 Cosmological mock data 

In order to assess the sensitivity of future cosmological observations to the axion mass we use 
the parameter forecast code 1 first described in reference [39] for galaxy clustering and cosmic 
shear surveys, and later extended in reference [38] to include cluster measurements. We adopt 
the same assumptions about the observational data, which conform to the specifications of 
the upcoming Euclid mission [37]. A summary of the mock data sets is as follows. 

1. The cosmic shear auto-spectrum C|® ■ runs in the range 2 < i < 7^ a x = 2000, and the 
indices i. j E [1,IV S ] label the redshift bin. As discussed in reference [39], shot-noise 
always dominates well before l reaches 2000, so that the precise value chosen for £ max 
is not critical. Redshift slicing is performed in such a way that all bins contain similar 
numbers of source galaxies and hence suffer the same amount of shot noise. We choose 
N s = 2, as no substantial improvement in the parameter sensitivities can be expected 
beyond this point [39]. 

2. The galaxy auto-spectrum Cf f ■ uses multipole moments running from l = 2 to fmLc i n 

redshift bins i,j E [l,IVg], where fflLx is set following the prescription of reference [39] 
so as to minimise the non-linear contamination. As shown in the said paper, most cos¬ 
mological parameter sensitivities saturate at around A r g = 11; we shall also adopt this 
value. The linear galaxy bias is assumed to be perfectly known. While this assump¬ 
tion is not very realistic, in the absence of a concrete bias model the only alternative 
would be to adopt the opposite extreme in which the bias is completely unknown and 
uncorrelated between redshift bins. We have tested this equally unrealistic assumption 
in a previous work [39], and found that the resulting parameter constraints are utterly 
uncompetitive if the galaxy data were used on their own, and, if analysed in combina¬ 
tion with shear data, hardly better than the shear-only constraints. We therefore do 
not consider the case of an unknown galaxy bias in the present work. 

3. The shear-galaxy cross-spectrum . in the shear redshift bin i E [1,IV S ] and galaxy 

redshift bin j E [1, N s ] runs from t = 2 to where £mi.x is determined by the galaxy 
redshift binning. 

4. The cluster mass function measurements IVjjVo source from clusters in the redshift 
range z E [0.01, Zhigh] and with masses M E [Mthr(z), 10 16 M 0 ], where Mthi(z) is the 
redshift-dependent mass detection threshold, and Zhigh the redshift at which 
exceeds 10 16 M@. The index i E [1, N z ] labels the redshift bins, with boundaries chosen 
such that for the fiducial model all bins contain equal numbers of clusters. In each 
redshift bin i we further subdivide the cluster sample into mass bins j(i ) E [l,IV m ], 
again demanding that the same number of clusters should fall into each bin [38]. 

5. A mock data set from a PLANCK-like CMB measurement is generated according to the 
procedure of reference [62], together with modifications according to reference [63] so 
that it better captures some of the gross features of the realistic Planck likelihood, such 
as the fact that its effective sky coverage is scale-dependent. 

1 http://jhamann.web.cern.ch/jhamann/simdata/simdata.tar.gz 
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3.2 Cosmological model space 

We use an 8-parameter cosmological model ©, consisting of six “vanilla” parameters ©( 6 ), 
extended by the sum of neutrino masses here denoted as m v and the axion mass m a , he., 

0 = = ^(wcdmj w b , h, ln(10 10 A s ), n s ),m v , m^J . (3.1) 

Here, w C( j m = H c d m ^ 2 and Wb = Hbh 2 are the present-day physical CDM and baryon densities 
respectively, h = Hq /(100 km/s/Mpc) is the dimensionless Hubble parameter, ln(10 10 A s ) and 
n s denote respectively the amplitude and spectral index of the initial scalar fluctuations, and 
z re is the reionisation redshift. 

For the part of the parameter space not related to neutrinos or axions, our fiducial 
model is defined by the parameter values 

©J2 = (0.0226,0.1126587, 0.7, 2.1 x 10“ 9 , 0.96, ll) , (3.2) 

while for m v and m a we shall test a variety of fiducial values. Specifically, for the axion mass 
we use a range of values 

m a ,m = 0.01,0.1,0.15,0.2, 0.5 eV, (3.3) 

covering the possibility of axion freeze-out both before and after the QCD phase transition 
(see figure 2). Neutrinos provide another source of hot dark matter of an unknown size. We 
expect their phenomenology to be similar to that of axions so that m„ and m a are likely to 
be correlated. Therefore, we vary also the fiducial neutrino mass sum in the range 

m u> m = 0.06,0.11,0.3, 0.7,1 eV, (3.4) 

i.e., from the minimal mass allowed in the normal hierarchy by neutrino oscillation experi¬ 
ments, to well beyond the current cosmological upper bound. 

Our analysis assumes, for the most part, a pessimistic scenario in which the neutrino 
mass sum is an unknown that has to be fitted simultaneously with the axion mass to the 
mock data, possibly leading to a degradation of the sensitivity to the latter. However, we 
shall consider also an optimistic case in which the neutrino mass is taken to be infinitely well 
measured by other means (e.g., by tritium /3-decay experiments); in practice this means we 
hold the neutrino mass sum fixed at its fiducial value when fitting the mock data. 

Lastly, we note that a number of recent cosmological analyses have found hints for a 
small hot dark matter contribution due to either a nonzero neutrino mass sum or axion 
mass [35, 64], These results are largely a consequence of tension between different observa¬ 
tions within the minimal ACDM interpretation—notably, the conflicting values of as and 
reported by cluster catalogues and by Planck’s CMB temperature anisotropy measurements— 
which tends to settle on a nonzero hot dark matter component as a middle ground. However, 
the tension between these observations may just as well indicate an incomplete understanding 
of the systematics [65]. 

Another example was the apparent conflict between the primordial tensor amplitudes 
inferred from Planck’s CMB temperature anisotropies and from the BICEP2 H-polarisation 
measurement. In this case, an increase to the universe’s dark radiation content [66, 67] had 
been proposed as a solution. However, in the end the discrepancy was shown to be due to 
insufficient modelling of the galactic polarised dust emission [52]. 

While the study of such tensions in actual data is interesting, for the purpose of this 
forecast we prefer not to artificially insert tensions into our mock data; we shall rather work 
with unbiased simulated observables, assuming that systematic issues in the measurements 
themselves will have been sorted out by the time the Euclid mission produces results. 
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all 


m„,fid = 0.06 eV 

to„, fid = 0.11 eV 

Ad = 0.3 eV 

TO^fid = 0.7 eV 

TO„, fid = 1 eV 

m a , fid = 0.01 eV 

0 - 0.136 

0 - 0.136 

0 - 0.144 

0 - 0.137 

0 - 0.144 

nVM = 0.1 eV 

0 - 0.144 

0 - 0.144 

0 - 0.144 

0 - 0.137 

0 - 0.144 

m a ,&d = 0.15 eV 

0.145 - 0.256 

0.145 - 0.281 

0.145 - 0.286 

0.145 - 0.303 

0.145 - 0.317 

to q ,fid = 0.2 eV 

0.155 - 0.302 

0.155 - 0.318 

0.158 - 0.352 

0.154 - 0.331 

0.155 - 0.359 

iria, fid = 0.5 eV 

0.415 - 0.570 

0.400 - 0.567 

0.383 - 0.566 

0.392 - 0.569 

0.341 - 0.581 

ccl 


ui„,fid = 0.06 eV 

to„, fid = 0.11 eV 

TO„, fl d = 0.3 eV 

TO^fid = 0.7 eV 

TO„, fid = 1 eV 

m aM = 0.01 eV 

0 - 0.144 

0 - 0.144 

0 - 0.144 

0 - 0.135 

0 - 0.144 

m a fi d = 0.1 eV 

0 - 0.136 

0 - 0.138 

0 - 0.137 

0 - 0.137 

0 - 0.137 

m a , fid = 0.15 eV 

0.145 - 0.392 

0.145 - 0.507 

0.145 - 0.749 

0.145 - 0.706 

0.145 - 0.893 

uia,fld = 0.2 eV 

0.157 - 0.407 

0.161 - 0.531 

0.168 - 0.825 

0.169 - 0.869 

0.173 - 0.902 

m aM = 0.5 eV 

0.233 - 0.642 

0.204 - 0.765 

0.266 - 0.933 

0.226 - 1.033 

0.253 - 1.107 

csgx 


m vM = 0-06 eV 

m V) fld = 0.11 eV 

TO^fid = 0.3 eV 

TO„, fld = 0.7 eV 

TO„, fid = 1 eV 

m a , fid = 0.01 eV 

0 - 0.137 

0 - 0.137 

0 - 0.144 

0 - 0.136 

0 - 0.137 

m aM = 0.1 eV 

0 - 0.144 

0 - 0.136 

0 - 0.144 

0 - 0.137 

0 - 0.137 

m a ,&d = 0.15 eV 

0.145 - 0.309 

0.145 - 0.309 

0.145 - 0.355 

0.145 - 0.331 

0.145 - 0.354 

m aM = 0.2 eV 

0.151 - 0.311 

0.155 - 0.330 

0.153 - 0.387 

0.154 - 0.352036 

0.157 - 0.381 

m a ,&d = 0.5 eV 

0.412 - 0.570 

0.396 - 0.567 

0.361 - 0.566 

0.375 - 0.563 

0.278 - 0.584 

CS 


uiy.fld = 0.06 eV 

fid = 0.11 eV 

TO„, fid = 0.3 eV 

TO„,fi d = 0.7 eV 

TO^d = 1 eV 

m a fid = 0.01 eV 

0 - 0.137 

0 - 0.137 

0 - 0.144 

0 - 0.144 

0 - 0.144 

m a , fid = 0.1 eV 

0 - 0.144 

0 - 0.137 

0 - 0.144 

0 - 0.137 

0 - 0.144 

m n ,sd = 0.15 eV 

0.145 - 0.321 

0.145 - 0.350 

0.145 - 0.513 

0.145 - 0.704 

0.145 - 0.824 

m a , fid = 0.2 eV 

0.152 - 0.349 

0.157 - 0.374 

0.161 - 0.528 

0.164 - 0.738 

0.166 - 0.838 

m a ,m = 0.5 eV 

0.200 - 0.585 

0.186 - 0.604 

0.208 - 0.711 

0.187 - 0.865 

0.218 - 0.941 


Table 1. ID marginal 95% credible intervals for the axion mass inferred from four different data 
combinations (indicated above the panels; see text for an explanation of the abbreviations), assuming 
different fiducial m a and m„ values. 


4 Numerical results 

4.1 Unknown neutrino mass 

Table 1 shows the allowed ranges of axion masses for our 25 different combinations of fiducial 
axion and neutrino mass values, inferred from four combinations of data set: 

• ccl: CMB + clusters, 

• cs: CMB + shear auto-spectrum, 

• csgx: CMB + shear auto-correlation + galaxy auto-correlation + shear-galaxy cross¬ 
correlation, and 
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Figure 3. Left: Scatter plot of inferred mass values in the m a -m v plane, assuming a common 
to i/.fid = 0.7 eV, and m a , fid = 0.15 eV > toqcd (top) and m a ,fid = 0.10 eV < tttqcd (bottom). 
The combination of all data sets is considered in both cases. The colours of the points indicate 
the corresponding axion number density according to the colour bar. Right: ID marginal posterior 
probability density distribution in m a inferred from various combinations of data sets, assuming 
m^fid = 0.06 eV and two choices of fiducial axion masses below the QCD mass limit. 


• all: CMB + shear auto-correlation + galaxy auto-correlation + shear-galaxy cross 
correlation + clusters. 

The “jump effect” discussed in section 2.2 is immediately apparent: for fiducial axion 
masses below ttiqcd = 0.145 eV, it is not possible to diagnose a nonzero axion mass even with 
a future mission such as Euclid. This is in sharp contrast to the case of neutrino masses, 
where it has been shown that even the minimum mass sum of m u = 0.06 eV can be detected by 
Euclid with high statistical significance [39]. As explained earlier, this is because axions with 
masses smaller than ttoqcd decouple at temperatures above Tqcd = 170 MeV, so that their 
corresponding present-day number densities suffer from entropy dilution acquired during the 
QCD phase transition, as already illustrated in figure 2; in contrast, the present-day neutrino 
number density is always n u = 112 cm= 3 per flavour, independently of the neutrino mass. As 
soon as the fiducial axion mass exceeds mQCD, however, a positive detection becomes easily 
possible. 

We further illustrate this jump effect in the left panel of figure 3, which shows two scatter 
plots of the m a and m v values from the likelihood analysis of the “all” data combination, 
assuming a common m Ut fid = 0.7 eV but two different m a fi d values: one just above tuqcd in 
the top panel (m Qi gd = 0.15 eV), and one just below in the bottom {rn a fi& = 0.1 eV). The 
axion number density corresponding to each scatter point is denoted by its colour coding. 
Clearly, in the top panel, the likelihood analysis returns a high concentration of points at the 
fiducial axion mass value. In the usual fashion the rest of the points are scattered around 
m a , fid, but with the twist that the density of points is sharply reduced for axion masses 
below niQCDj the QCD phase transition mass limit. 
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Figure 4. 2D marginal 68% and 95% credible contours for m a and m u derived from the “ccl” (red), 
“cs” (green), and “all” (black) data combinations. The matrix of fiducial mass combinations is the 
same as in table 1. 

A similar cut-off can also be seen in the bottom panel, except that in this case, because 
m a< fid < U7.QCD) it is those axion masses exceeding ttiqcd that are devoid of scatter points, all 
of which are now spread instead between 0 eV and tuqcd- Indeed, the probability distribution 
between m a = 0 eV and tu.qcd is essentially flat in the m a -direction no matter the exact 
choice of m a , fid; as demonstrated in the right panel of figure 3 by the ID marginal posteriors 
for a selection of m a ^ values. This flatness signifies that the m a% m < "Tqcd cases are for 
all purposes indistinguishable from one another. 

Lastly, figure 4 shows the 68% and 95% contours in the 2D (m a , ?n I/ )-parameter space 
derived from various data combinations and for the full matrix of fiducial mass values shown 
in table 1. For most fiducial mass combinations there is a strong degeneracy between 
and m a when the “ccl” and “cs” data combinations are considered separately. This degen¬ 
eracy is however almost entirely broken for the “all” data combination, where the galaxy 
auto-spectrum is also fitted along with shear and cluster data. To understand why this hap¬ 
pens, consider the scatter plot in the 3D (Hq, m a , m^-parameter space displayed in figure 5. 
Evidently, the degeneracy directions resulting from the “cs” and “cs+clusters” data combi¬ 
nations are close to orthogonal to that from the galaxy auto-spectrum alone. Thus, adding 
the galaxy auto-spectrum not only breaks the degeneracy between m a and m u , but also the 
two parameters’ respective degeneracies with Hq. 
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Figure 5. Scatter plot in the 3D (Ho, m a , TO„)-parameter space for the fiducial mass combination 
m a ,fid = 0.2 eV and d = 0.3 eV. The orange points denote the “cs” data combination, the green 
points “cs+clusters”, and the black points the galaxy auto-spectrum. The MCMC fitting has been 
performed only in this 3D subspace, while all other parameters have been held fixed at their fiducial 
values. For improved readability we plot only those points satisfying Ay 2 < 2. 

4.2 Fixed neutrino mass 

Cosmological data of the accuracy studied here will not be available until about a decade 
from now. It is therefore conceivable that the neutrino mass will have been pinned down 
before then, via, e.g., measurements of tritium decay by the KATRIN experiment or by a 
post-KATRIN generation experiment [68] . We therefore study also the extreme case in which 
the neutrino mass is assumed to be exactly known, in order to assess the maximum possible 
effect our prior knowledge of neutrino masses has on the cosmological sensitivity to axion 
masses. We consider the same 25 combinations of fiducial axion and neutrino masses, but 
vary now only the axion mass as a fit parameter. The likelihood analysis is restricted to the 
“all” combination of data sets. 

Table 2 shows the allowed m a ranges for all 25 combinations of m a ,m and m u $ d- Com¬ 
paring with their counterparts in table 1, we can immediately make the following observations: 
On the one hand, for fiducial axion masses below the QCD mass limit, what happens to the 
neutrino mass continues to have virtually no impact on the sensitivity to m a . On the other 
hand, in those cases where m a ,fid exceeds jtiqcd, the corresponding error bands shrink by 
~ 30% to ~ 60% upon fixing the neutrino mass. 

These results can be understood from an inspection of the 2D marginal credible contours 
in figure 4. When the neutrino mass is left as a free parameter, construction of the ID 
marginal posterior for the axion mass m a consists essentially of integrating the 2D ellipse 
along the jn^-direction. With the neutrino mass fixed, however, the ID marginal posterior 
for m a is simply the cross-section of the ellipse at m v = d- Thus, because the ellipses 
in the m a ,fid < uiqcd cases are almost exactly aligned with the m a -axis, the two procedures 
yield very similar results. In contrast, the ellipses in the m a ^ > mQCD scenarios are all 
inclined at an angle; assuming that the ellipses represent perfect 2D Gaussian distributions, 
integration along the m^-direction will always yield a ID posterior for m a that is broader than 
the cross-section at m u = rri v ^ ( \, with the greatest disparity occurring when the inclination 
approaches 45°. 
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all, fixed neutrino masses 


m Vj dd = 0-06 eV 

m v , fid = 0.11 eV 

m v fi d = 0.3 eV 

m v> fid = 0.7 eV 

m Ui fl d = 1 eV 

m a , ad = 0.01 eV 

0 - 0.135 

0 - 0.134 

0 - 0.136 

0 - 0.137 

0 - 0.137 

»Vfid = 0.1 eV 

0 - 0.138 

0 - 0.137 

0 - 0.137 

0 - 0.144 

0 - 0.137 

m-afid = 0.15 eV 

0.140 - 0.211 

0.140 - 0.213 

0.140 - 0.212 

0.140 - 0.215 

0.140 - 0.223 

uio.fid = 0.2 eV 

0.159 - 0.250 

0.158 - 0.253 

0.157 - 0.253 

0.157 - 0.256 

0.154 - 0.256 

m a ,M = 0.5 eV 

0.454 - 0.544 

0.453 - 0.545 

0.454 - 0.544 

0.455 - 0.543 

0.454 - 0.543 


Table 2. ID marginal 95% credible intervals for the axion mass inferred from the “all” data combi¬ 
nation, assuming different fiducial m a and m v values. In each case the neutrino mass has been fixed 
at the fiducial value in the analysis. 


4.3 When can a detection be claimed? 

Having calculated the formal sensitivity to m a for a variety of different fiducial models we 
might now also ask the question: How strongly would a model with m a = 0 be disfavoured 
depending on the assumed fiducial model? 

Table 3 shows the Ay 2 values for the best-fit model assuming m a = 0, relative to the 
global best fit (i.e., the fiducial model), derived from the “all” data set for the 25 fiducial mass 
combinations. (Note that we have also varied the neutrino mass as a fit parameter.) Looking 
at the first two rows, it is immediately clear that fiducial axion masses below 0.15 eV, for 
which Ay 2 ~ 0, will not in practice be distinguishable from zero given the accuracy of the 
data used here. This conclusion is consistent with our understanding of the Bayesian results 
presented in section 4.1 and 4.2. 

As soon as the fiducial axion mass surpasses 0.15 eV, the A% 2 value jumps to the 
range 13-37, depending on the fiducial neutrino mass; the larger the value, the lower 

the A\- Physically, this trend can be traced to a very small residual degeneracy between 
m a and m u , so that a 1 eV neutrino could partially fill in the role of a 0.15 eV axion while 
a 0.06 eV neutrino could not. But this is ultimately only an academic point of interest, 
because A\ 2 values as large as 13-37 at the cost of just one additional fit parameter would in 
any case constitute “strong” evidence in favour of a nonzero axion mass, no matter whether 
we apply the Bayes information criterion (BIC), the Akaike information criterion (AIC), or 
other constructs. 

We therefore conclude that for any axion mass above the mass threshold 0.15 eV re¬ 
quired for the axion to decouple after the QCD epoch, Planck data in combination with a 
EuCLiD-like photometric survey will have sufficient accuracy to reliably measure the axion 
mass. If on the other hand the real axion mass falls below the QCD mass threshold, then 
the corresponding axion number density is so low that a future survey with a much better 
accuracy would be needed to further tighten the axion mass bound from above. 

5 Conclusions 

The next generation of large-volume photometric surveys, notably Euclid, together with 
CMB data from the Planck mission are almost guaranteed to reach a sensitivity to the cosmic 
hot dark matter content good enough to measure the sum of neutrino masses m v with 
high significance. This is so even if m u should take on its most pessimistic value, 60 meV, 
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m uM = 0.06 eV 

m l , i ad = 0.11 eV 

m V} ad = 0.3 eV 

m Vi M = 0.7 eV 

mvM = 1 eV 

mafid = 0.01 eV 

0 

0 

0 

0 

0 

m a , fid = 0.1 eV 

0 

0 

0 

0 

0 

m a ,fid = 0.15 eV 

37 

36 

31 

18 

13 

ni a fid = 0.2 eV 

47 

46 

40 

25 

18 

m a ,M = 0.5 eV 

75 

71 

64 

55 

38 


Table 3. Ay 2 values (rounded to the nearest whole number) for the best-fit model assuming m a 
fixed at zero, relative to the global best-fit, derived from the “all” data set for the 25 fiducial mass 
combinations. 


the minimum established by flavour oscillation experiments. Complementarily, laboratory 
experiments will provide a direct measurement of the effective electron neutrino mass, albeit 
likely on a longer timescale. Depending on what is actually found at these experiments, 
there will be a new frontier of either identifying an additional hot dark matter component, 
or placing very precise constraints on non-neutrino contributions, e.g., in the form of a hot 
axion population. 

Motivated by the IAXO proposal for a next-generation solar axion search, we have 
explored the range of axion masses that is accessible to cosmological hot dark matter searches. 
For large enough axion masses, thermal axions are produced after the QCD epoch by very 
efficient axion-pion interactions, so that the resulting axion population is almost comparable 
to that of one species of neutrinos. For small masses, however, axion freeze-out occurs before 
the QCD epoch, where axion production from interactions with free quarks and gluons is 
much less efficient, and the final axion population is invariably diluted by the large number 
of colour degrees of freedom that disappear when the temperature of the universe drops 
below Tqcd- 

Based on this understanding, and assuming that Tqcd = 170 MeV, we find that axion 
masses larger than ~ 0.15 eV can be easily pinned down by cosmology. If on the other hand 
the axion mass is lower than this threshold, then cosmological observations are essentially 
blind to the hot axion population for the foreseeable future. 

When optimising the IAXO apparatus and search strategies, we conclude that the mass 
range rn a < 0.15 eV should have the highest priority because larger axion masses can be 
constrained by cosmology. To study this higher mass range with IAXO requires helium filling 
of the conversion region to achieve a match between the axion search mass and the refractive 
photon mass, a method that has been used at CAST and requires sequential pressure settings 
to cover a broad range of axion search masses. Therefore, cosmology can help to simplify the 
IAXO search strategy in that the technically most difficult mass range may well have been 
excluded by cosmology before IAXO becomes operational. 
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